Correlations in optically-controlled quantum emitters 
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We address the problem of optically controlling and quantifying the dissipative dynamics of quan- 
tum and classical correlations in a set-up of individual quantum emitters under external laser ex- 
citation. We show that both types of correlations, the former measured by the quantum discord, 
are present in the system's evolution even though the emitters may exhibit an early stage disen- 
tanglement. In the absence of external laser pumping, we demonstrate analytically, for a set of 
suitable initial states, that there is an entropy bound for which quantum discord and entanglement 
of the emitters are always greater than classical correlations, thus disproving an early conjecture 
that classical correlations are greater than quantum correlations. Furthermore, we show that quan- 
tum correlations can also be greater than classical correlations when the system is driven by a laser 
field. For scenarios where the emitters' quantum correlations are below their classical counterparts, 
an optimization of the evolution of the quantum correlations can be carried out by appropriately 
tailoring the amplitude of the laser field and the emitters' dipole-dipole interaction. We stress the 
importance of using the entanglement of formation, rather than the concurrence, as the entangle- 
ment measure, since the latter can grow beyond the total correlations and thus give incorrect results 
on the actual system's degree of entanglement. 

PACS numbers: 03.65.Ud, 03.67.Mn, 42.50.Fx, 03.65.Ta, 34.80.Pa 



I. INTRODUCTION 

The nature and the optimal dynamical control of the 
correlations present in a given interacting quantum sys- 
tem are of interest to fundamental aspects of quantum 
physics but also to the physical implementation of quan- 
tum computation and algorithms A decade ago, 
two separate works reported from Los Alamos [3] and Ox- 
ford [5[ set a formal framework that allowed the identifi- 
cation of quantum correlations different to entanglement. 
Such correlations ubiquituosly manifest themselves in in- 
teracting quantum systems, even in the absence of en- 
tanglement, and it has been demonstrated that they can 
play a key role in quantum computing protocols such as 
DQC1 [y, [7|, without the need to invoke quantum en- 
tanglement. It has become clear that zero entanglement 
does not imply classicality, and the task of classifying 
the kind of total correlations present in a given system 
into quantum and classical parts has led to the introduc- 
tion of quantifiers (different to entanglement) such as the 
quantum discord (QD) Q. This measure of the 'quan- 
tumness' of correlations is defined in terms of the differ- 
ence between the two possible ways to obtain the mutual 
information (MI) that arises in the quantum case. 

More recently, it has been argued that the power of 
the quantum computer can be attributed to both en- 
tanglement and quantum discord which are distributed 
following a monogamic relation 0. Further, other re- 
ports on the classification and quantification of quantum 
and classical correlations have been put forward [j^-fli]). 
A unified view of quantum, classical and total correla- 



tions using the concept of relative entropy as a distance 
measure of correlations has been discussed in [Til- An 
experimental investigation of classical and quantum cor- 
relations has been recently reported for a system of two 
polarized photons in the presence of a phase-decoherence 
non-dissipative channel |l5| , and a sudden transition be- 
tween classical and quantum correlations dynamics has 
been reporte d 1151 , a result in line with a recent theoret- 
ical proposal [12j. 
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We start with a brief definition of correlations (Sec. HI]) 
followed by the formulation of the dissipative dynamics 
under consideration, in Sec. IIIII Our main findings are 
detailed in Sections IIVIIVII as follows: We show that the 
entanglement of formation, rather than the concurrence, 
is the free of ambiguity metric that should be used to cor- 
rectly quantify the entanglement generated in the emit- 
ters' system (Sec. IIV|) . The computed concurrence can 
reach higher values than the mutual information, which 
is the highest possible value for the total (classical plus 
quantum) correlations. In Sec. |V] we demonstrate ana- 
lytically, in the absence of laser pumping, that there exist 
an entropy bound for which the quantum correlations are 
always greater than classical correlations, a fact which 
disproves an early conjecture posed in [25|. In Section 
IVII we give numerical examples for which such result also 
holds when the emitters arc optically-controlled by an 
external laser field. This also gives a method for station- 
ary enhancement of the degree of correlations by tailoring 
the laser amplitudes and the strength of the dipole-dipole 
interaction. 
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II. QUANTIFYING CLASSICAL AND 
QUANTUM CORRELATIONS 

We begin by considering an open bipartite system, say 
the molecular dimer AB. Its quantum dynamics can be 
described in terms of the total density matrix pab] the 
partial information about each of the subsystems is given 
by the reduced density matrices pi = trj(pAB), hi = 
A,B. 



A. Mutual Information 

In classical information theory, the mutual informa- 
tion (MI) can be written in terms of the conditional 
entropy of two random variables, that is, for a system 
described by two classical random variables, the most 
information that one can gain from the two variables 
is the MI. Let X and Y be two random variables with 
Shannon entropies H(X) = — J2 x P(x=x)^ogP(x=x), and 
H(Y) = — J2 y P(Y=y)\°EP(Y=y)- The MI is written as: 

I(X :Y)=H(X) + H{Y)-H(X,Y). (1) 

Using Bayes's rules for the conditional probability 

P(x\Y=y) = P(x,Y= v )/P(Y=y), the MI reads 

J(X:Y)=H(X)-H(X\Y), (2) 
where H(X\Y) is the conditional entropy (l6j . 

B. Total correlations 

The mutual information has been demonstrated to de- 
scribe the whole content of correlations in a given quan- 
tum system, and at any time t, the total correlations 
present in such a system can be quantified by the quan- 
tum version of Eq. (fTJ), the so-called quantum mutual 
information 

I(pab) = S(p A ) + S(p B ) - S(pab), (3) 

where S(p) = —tr(plog 2 p) is the von Neumann entropy. 
The problem that arises in the quantum case is that the 
non trivial version of Eq. is not always identical to 
I(pab)- This is because, quantum mechanically, the con- 
ditional entropy implies that the state of a subsystem 
(say A) is only known after a set of measures has been 
carried out on the other subsystem (B). Ollivier and 
Zurek Q proposed a generalization of Eq. © assum- 
ing a complete unidimensional projective measurement 
made on subsystem B, with projectors {n^} satisfying 

E(nf ) t nj 5 = 1. The generalization of Eq. © allows 
another way of writing the quantum mutual information 
(which in general does not coincide with that of Eq. © ) : 

J(PAB){af} = S(pa) - S(p Al{Tl B } ), (4) 



where the conditional entropy S(pa\b) = S(Pa\{yi b }) = 
12jPj S (pA\n? )> with probability p 3 = tr(Uf p A B^lf), 
and the density matrix after the measurements have been 
performed on subsystem B is given by 

p ^f ~ HTifpABiify (5) 

Equation © gives the amount of information gained 
about A after a measure has been carried out on B. 



C. Classical correlations 

A measure for classical correlations (CCs) was intro- 
duced in Ref. This quantity is defined as the maxi- 
mum extractable classical information from a subsystem, 
say A, when a set of positive operator valued measures 
(POVMs) [lH has been performed on the other subsys- 
tem (B): 

CC(pab) = max jn B) J(p A B){nf} = ( 6 ) 

max {ri B } (s{p A ) - ^PjS(p A )^j , 

j 

where S(p 3 A ) is the entropy associated to the density ma- 
trix of subsystem A after the measure. Briefly, a measure 
of classical correlations should incorporate the following: 

i) CC(pab) = if and only if pab = Pa <8> Pb, and 

ii) CC must be non-increasing, and invariant under lo- 
cal unitary operations. In addition, the set of POVMs 
that maximize the classical correlation Eq. ([5]) is a com- 
plete unidimensional projective measurement jTl^} [13]. 
If {|0), |1)} defines the basis states for qubit B, the pro- 
jectors can be written as 11^ = 1 <E> \j) j = a, b, where 
| a) = cos(9|0) + e^sin(9|l), |6) = e-^sin0|O) - cos0|l), 
and the optimization is carried out over angles and <fi. 

D. Quantum correlations 

1. Quantum Discord 

Following Sec. Ill 131 and given that Eqs. © and (gj 
are not always equivalent, a measure of the quantumness 
of correlations — the so-called quantum discord (QD) is 
defined as the minimum of the difference between the 
Eqs. © and ©. The QD then reads 

D(pab) = min {n s } \ I(pab) - J(pab){ U ? }) = ( 7 ) 

S(p B ) - S(pab) + min { nf } ^2PjS(p A \nf)> 

j 

and the measure of total correlations becomes I (pab) = 
D(pab) + CC(pab), where D(pab) is a measure of 
'purely quantum' correlations, and CC(pab) is the max- 
imum extractable classical information. Briefly, some 
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properties can be inferred from Eq. ([7]): i) For pure 
states, D = S(pb), ii) £) = if and only if pab = 
V, ll^/M/;llj ; - iii) < D < S{p B ), iv) D is an entan- 
glement monotone on pure states. Both QD and CC are 
antisymmetric by definition, i.e., they depend on what 
subsystem the measures are taken on [1 01 III. With- 
out loss of generality, we take qubit B to be the measured 
one. 



2. Quantum Entanglement 

Entanglement is a special and a very particular kind of 
correlation [l| , radically different to the quantum discord. 
It is fundamental in quantum phenomena, and aprecious 
resource for performing quantum computing [2 Q and 
quantum communication protocols (l6j . A bipartite pure 
state is said to be entangled if at least two cofficients of its 
Schmidt decomposition do not vanish [l6| , and the opera- 
tors pa and ps have the same non- vanishing eigenvalues: 
they are equal to the square of the Schmidt numbers. A 
state acting on Hilbcrt space T-Lab is separable if it is of 
the form p = J2i=i c iPt ® pf> f° r some k, where pf and 
pf are states on %a and T-L b respectively. If p is a pure 
state, i.e., p = \&ab) {^ab\, then it is separable if and 
only if \^ab) = \^a) ® |*b). 

The entanglement of formation (EoF), arguably the 
most prominent measure of entanglement, allows an an- 
alytical expression in the bipartite case. The EoF for the 
density matrix p is defined as the average of the entan- 
glement of all the states composing p, minimizing over 
the whole composition (over pure state decompositions) 

m, 

EoF(p)=minJ2Pi^i), (8) 

i 

where E(ip) = -tr[p A log 2 pa] = -tr[p B log 2 Pb] is the 
entropy of each of the subsystems. The EoF can also 
be written in terms of the quantum mutual information 
as EoF {pab) = min^pi \l{p\ '■ Pb)- Although the 
minimization is not trivial in general, the EoF has an 
analytical solution for the two-qubit case [HI ■ 
For pure states, the EoF can be written as 

EoF{i>)=£{C{4>)), (9) 

where the concurrence C{ip) =| (ip\ip) |, and the function 

where h{x) — —x log 2 x — (1 — x ) log 2 (l — x) denotes the 
binary entropy function. £ is monotonically increasing, 
and it goes from to 1 when the concurrence C ranges 
from to 1. This is why C can be used as an entangle- 
ment quantifier. This said, it should be stressed that it 
is only the EoF that is an entanglement measure, and C 
gets its meaning via its relation to the EoF. 



For mixed states, the EoF is written as [l8| 

EoF{p)=£{C{p)). (11) 

Consider, in decreasing order, the eigenvalues of the 
matrix V Pab Pab, where pab = {o y (8 o- y )p A B{o- y ® cr v ), 
Pab is the elementwise complex conjugate of p, and o~ y 
is the Pauli matrix. The concurrence C is defined as 

C{p AB ) =max{0,A 1 - A 2 - A 3 - A 4 }, (12) 

where the A^'s are as defined above or, equivalently (also 
in decreasing order), the square root of the eigenvalues 
of the non-Hermitian matrix pab Pab (HI- 

III. SINGLE QUANTUM EMITTERS 

We consider dipole-dipole (d-d) interacting quantum 
emitters in contact with a surrounding medium that 
acts as an environment. This is the case, for exam- 
ple, of coupled terrylene single molecules embedded in 
a para-tcrphenyl crystal for which symmetric and anti- 
symmetric (super- and sub-radiant) states, as well as d-d 
interaction strength, and single and collective decay rates 
have been measured [l||. For these, molecular quantum 
control and entangling quantum logic gates have been 
designed [20l |. in addition to the control of the quantum 
dissipative entanglement dynamics [U, [24|. Moreover, 
quantum control experiments and femtosecond single- 
qubit gates involving single organic molecules at room 
temperature have recently been performed [22[ ■ 

For the purpose of the analysis of the dynamics of the 
correlations presented here, we consider that the emitters 
couple with interaction strength V, and have individ- 
ual transition frequencies fj (molecule i) in interaction 
with the quantized radiation field. We denote by |0j), 
and |1;), i = 1,2, the ground and excited state of emit- 
ter i, respectively, and hence the single particle Hamil- 
tonian Hq = —^vicr^ — ^2°^ ■ The two two- level 
emitters are separated by the vector ri 2 and are char- 
acterized by transition dipolc moments fii = (Oj|Dj|l,-), 
with dipole operators , and spontaneous emission rates 
r^. The emitters arc externally controlled by a coherent 
driving laser field which acts on each of the molecules 
with a coupling amplitude Mi = —fMi-Ei, and a fre- 
quency u>l- The light-matter interaction Hamiltonian 
H L = h£W(a { i ) e i,JJLt + o- ( _i ) e- iu ' Lt ), where Ei is the ampli- 
tude of the coherent driving acting on molecule i located 
at position r,-. We also allow for the molecular and laser 
detunings A_ = Ui — f 2 , and A + = Ul j^ 2 — v^. 

IV. DISSIPATIVE DYNAMICS OF QUANTUM 
CORRELATIONS 

The emitters interaction Hamiltonian can be written in 
the computational basis of direct product states \i) ® \j) 
{i,j = 0,1) as Hs = H o + ^12, where H12 is set 
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by the dipole coupled molecules of interaction energy 



tt _ hV((l) 
tl\2 — —\<Jx 



< 2 > 



Hence, the dis- 



sipative dynamics corresponding to the ligth-matter in- 
teraction [23l |24j can be described, with H = Hs + ^l, 
by the quantum master equation 



P = -j j [H,p]+L{p), (13) 
where the dissipative Lindblad super-operator reads [2l[ 



L{p) = ~f {pa^^ 



2 ( xf 1) M 1) 



(2) (2) „ „ (2) - (2) 

u\'u_ p — 2a pa\' 



12 (1) (2) 
po\ a_ 



r 

~2 

T21 / . (2) (1) 



(i) (2) „ 

<T + <J_ p 



O (!) - (21 

2a_ pa\ 



(2) (1) - „ (2) „ (I] 
Cl <7_ p — 2(71 jOCTl 



where er^ = |lj)(0i|, and <7_ = |0j)(li| arc the rais- 
ing and lowering Pauli operators acting on emitter i. Ti, 
and Ti2 = = 7 are the individual, and the collec- 
tive spontaneous emission rates, respectively. The latter 
arises from the coupling between the emitters through 
the vacuum field. The interaction strength and the inco- 
herent decay rate read 



V = - [Ai ■ fi 2 - (Ai • ri2)(#2 • fi 2 ) 

[Ai ■ A2 - 3(Ai • ?i2)(A2 • fi 2 )] 

7 = c([Ai • A2 - (Ai • ria)(A2 • fi2)] 
[Ai ■ A2 - 3(Ai ■ ri 2 )(A2 • ri 2 ) 



cos z 



(14) 



sinz 

z 

cos z 



(15) 



smz 



where z = nk^r^, n denotes the medium refraction in- 
dex, k = ^,u - °™ 1 p - 3Vrrr * 



andC 



A. Identical emitters without laser field 

We begin with the simplest case scenario, that of iden- 
tical emitters or qubits, A_ = 0, no external laser field: 
ti = 0, and A + = [y x + v^)j2. 

We consider the emitters dipolar moments to be par- 
allel amongst them and perpendicular to the interqubit 
position vector 172 . In particular, if m Ao, the 
maximum (minimum) interaction strength and decay 
rate are obtained for parallel (perpendicular) dipole mo- 
ments as hV - 3h £y* [Ai ' A2 - 3(Ai • r u )(A2 ■ fi 2 )], 
7 = VT1T2 Ai ' A2- However, our calculations make di- 
rect use of Eqs. (Tl"4"|) and ([T5|) . since we investigate the 
dependence of correlations on the interqubit separation. 

The effect of dissipation on the correlations dynamics is 
initially shown in Fig. [T]for the decay rate 7 = 0.6884 V, 
and Ti = r 2 = T. Here, weakly coupled qubits are con- 
sidered for interaction strengths of the order of the decay 
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FIG. 1. (Color online) Correlations dynamics for a pair 
of interacting emitters initially in the doubly-excited state 
1 11). MI (solid thick line, purple), QD (solid thin line, blue), 
CC (doubly dashed line, green), C (dotted line, pink), EoF 
(dashed line, red). 7 = 0.6884 V. Main: T = 0.7818 V 
(ri2 = Ao/8). Inset: Functional relation between EoF and C. 
These functions are monotonically increasing from to 1. The 
zoom points out the differences between the EoF (G [0, 0.025]) 
and C (£ [0,0.1]). 



rates V ~ T which in the main graph corresponds to 
fi2 = Ao/8. The correlations last longer times and reach 
higher values the smaller T. As expected, the quantum 
mutual information is always greater than the rest of the 
correlations. For the emitters in the initial doubly excited 
state 1 11) we find that the entanglement quantifiers EoF 
and C exhibit a 'sudden birth' of entanglement [2l], [27j 
at a time r ~ 5/V (main graph); this behavior is 'short- 
ened' as the decay rate increases (not shown). In contrast 
to this, the quantum discord is always greater than en- 
tanglement (measured by EoF and C), and only vanishes 
around a time 3/V: only around a small time neighbour- 
hood centered at this value do the classical correlations 
(CC) become slightly greater than quantum correlations 
(QD). For the remaining parts of the evolution, QD > 
CC. We shall return to this point later, where we address 
the question of whether it is at all feasible to characterize 
a dynamical order relation between classical and quan- 
tum correlations. We remark that even in the absence 
of entanglement (f < r) there are both, quantum and 
classical correlations present in the dynamics, and that 
the QD is the only quantum correlation present in such 
a time frame. 



B. Entanglement of Formation vs Concurrence 

The concurrence has been widely used as an entan- 
glement metric due to the fact that, according to Eq. 
([TT]) . the function EoF(C) is monotonically increasing 
and ranges from to 1 as C goes from to 1. This said, 
the difference in value between these quantifiers can be 
significant and this is reflected in the interpretation of the 
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FIG. 2. (Color online) Same parameters and notations as in 
Fig- [Ij(main), but in the presence of an external continuous 
laser field of amplitude I = 0.4F in resonance with the qubits 
transition energies: A+ = (uil = wi = <^2 = w). The initial 
condition is taken to be the doubly-excited state |1I). 



correlation dynamics. The EoF changes at a slower rate 
than C at values close to 0; for instance, in the interval 
C € [0,0.1], the EoF is much smaller than C, as can be 
seen in the zoom of the inset of Fig. [1] The whole range 
for the metric is plotted in the inset of Fig. [TJ where 
the dashed line is a guide to the eye and indicates the 
equality EoF = C, to better appreciate the differences 
between EoF and C (solid curve). 

As for the physical interpretation of the correlations 
dynamics, we relate the above discussion with that pre- 
sented in Fig. [TJ Here, the dotted (pink) curve denot- 
ing C gives the impression that entanglement is almost 
the only quantum correlation present in the system for 
times around or greater than twice the sudden birth time 
(t > 2t), since for this, C and QD are almost, or actu- 
ally, superimposed. The solid (red) curve for the EoF 
clearly shows that most quantum correlations, however, 
in the system are far from entanglement. As a matter 
of fact, entanglement, as measured by EoF, contributes 
very little to the full quantum correlations dynamics. 

This point is made even more evident when we con- 
sider, for example, the influence of an external continu- 
ous wave laser field of amplitude I, and frequency ljl in 
resonance with the qubits transition energies: lul = u 
(below we shall return to the discussion of the laser ef- 
fects), as shown in Fig. [2J The physical parameters and 
the notations are the same as those shown in the main 
graph in Fig. [TJ The result that marks the importance of 
using EoF as the metric of entanglement is that, while in 
Fig. [2]the EoF (red curve, just below the discord which is 
plotted in blue) shows a behavior that can be understood 
following the lines of reasoning in Fig. [JJ (but now in the 
presence of a laser field; note the differences), the con- 
currence (dotted pink curve) exhibits a 'birth' of entan- 
glement (following the qualitative behavior of the EoF), 
just around a time that marks the beginning of the dom- 
inance of classical over quantum correlations, that grows 
well beyond the mutual information: since the latter is, 
by definition, the measure of total correlations present 



in the system at any given time, concurrence cannot be 
a good quantifier of entanglement in this context. This 
is an important point to bear in mind when performing 
quantum information protocols, where entanglement is 
required as a physical resource. As already pointed out, 
it is only the EoF that is an entanglement measure, and 
the significance of C is actually revealed through its re- 
lation to the EoF. Hereafter, we shall consider the EoF 
as the quantifier of the emitters' entanglement. 



V. HIERARCHY FOR THE EMITTERS 
CLASSICAL AND QUANTUM CORRELATIONS 

Regarding the conjecture that classical correlations are 
always greater than quantum correlations, and the re- 
cent experimental evidence that the opposite situation 
can also arise for a specific scenario in bipartite entan- 
gled photons subjected to a decoherent source (thus dis- 
proving the conjecture) [l5j . we ask whether it is at all 
feasible, and under what circumstances, to find and order 
hierarchy for the emitters' classical and quantum corre- 
lations. 

The following analysis of dynamical behavior for the 
correlations is given in terms of the inherent physical 
properties of the emitters, such as the collective damping 
(decay) rate 7, and the d-d interaction V. We demon- 
strate, in a simple manner, that the dimer system can 
be naturally dominated by quantum correlations at all 
times, that is, that QD or EoF are greater than classical 
correlations during the whole time evolution. For this, 
we solve the emitters' master equation (|13[) for identi- 
cal qubits {v\ = V2), with emitters only coupled through 
their dipolar interaction (there is no laser excitation). 
The calculation gives, for the non-trivial density matrix 
elements, the following result: 



1 

Poo, 00 — 2° 



-(7+r)t 



2e (7+r)t _ 2/(1 _ c 2 7t) 



cos ip 



B 2 ^ - 1 



Poi,oi — ^ c 



-1'/ 



e -7t(! + e 27t + 2/(1 _ e 2 ^) cos cp) 
(2 — 4a) cos 9 + 4/ sin ip sin 9 



Pio.io — ^° 



,-rt 



e-T'(l + e 2 ^ + 2/(1 - e 27t ) cos^) 
+(2 — 4a) cos — 4/ sin ip sin 9 



Poi.io — Pio.oi 
1 



4 6 



-rt 



e -7« (1 _ e 2 i" + 2/(1 + e 27 ") cos ip) 
2i(2/ sin </? cos9 + (1 - 2a) sin 9)1 , (16) 



where / = ^(1 - a)a, 9 = 2Vt, ip is a relative phase, 
and a defines the set of initial states \ct v ) = y^lOl) + 
e ltp ^/l — a|10), < a < 1, a choice motivated from states 
that can be 'naturally' prepared in the quantum emitters, 
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FIG. 3. (Color online) Hierarchy of correlations: D(pab) > 
CC(pab) /or a/Z 1 x = exp[-(r + 7)*]. The top (solid) 
curve is the quantity 2mins j . {5(p A | j n s j)}, and the bottom 

(dashed) curve corresponds to the entropy S(pab)- a = 1/2, 
<^ = 0; S( P a) = S( P b). 



e.g., separable states a = 0, 1, and symmetric (a = 1/2, 
and y> = 0) and anti-symmetric (a — 1/2, and ip = 7r) 
entangled states. For the class of a v initial states, whose 
time-evolution has been derived in Eq. (|16[) . the density 
matrix adopts the structure 



P(t) 



( Poo,oo 

Poi.oi Poi,io 

Poi.io Pio.io 

V 00 



(17) 



and thus the correlations can be analytically calculated. 
In particular, we are interested in finding a dynamical 
order hierarchy amongst classical and quantum correla- 
tions. Equations ^ and ([7]) allow to show that there is 
an entropy bound: 



D(pab) > CC(p AB ), 



(18) 



which reads 



S(p B ) + 2min Bj .{S'(p j4 | {rl f 3} )} - S(p A ) 



S(pab) > 0. 

(19) 

When the system begins in the maximally entangled sym- 
metric Bell state a = ^, <p = 0, the density matrix be- 
comes independent of the dipolar interaction V , and 7 is 
just an additional decay rate that makes the system decay 
with the new rate T + 7; furthermore, S(pa) — S(pb)- 

The independence of p with respect to V can be 
formally understood by recalling that the considered 
Bell states are part of a bigger family of states — pm, 
which have maximally mixed marginals, i.e., trA{pAi) = 
tr B {pAi) = 5/2x2, with I 2X 2 being the identity 2x2 
matrix; hence S(pa) = S(pb) = 1- Such states can be 
written as p M = \ (I + Yh=i hi<7i <8 [13, IM , where Oi 
are the Pauli matrices and the coefficients hi € 7Z satisfy 
constraints such that pu is a well defined density opera- 
tor. It can be shown that for the whole family of density 
matrices pm, the Markovian time evolution given by Eq. 
(fT5|) is completely independent of the interaction V. 



Returning to our case, quantum correlations are always 
greater than classical correlations. This can be seen in 
Fig. [3j where we have plotted the entropies condition 
that make D(pab) > CC(pab)- With the parametriza- 
tion x = exp[— (r + 7)t], in Fig. [3] we show that the 
quantum discord is greater than classical correlations for 
the whole dynamics. The equality is only obtained at the 
initial time (pure state) and at the 'final' (t — > 00) time; 
in the latter, the system tends to a separable and un- 
correlated state. In a similar manner, we can also show 
that the entanglement of formation EoF(pab) is always 
greater than the classical correlations. 

The same hierarchy of correlations can be proven by 
direct calculation. For the density matrix ([T7| . the con- 
ditional entropy, after the measurement, admits an ana- 
lytical expression and therefore both CC and D can be 
analytically computed: the conditional entropy has two 
possible solutions: 



S W (PA\{uf}) = -Ax),ook>g 2 
-Pio,iolog 2 



Poo, 00 



Poo, 00 + PlO.10 

PlO.10 
Poo, 00 + Pl0,10 



S {2) {PA\{Tlf }) = -5(1 - £)l0g 2 - 



~(l + f)Iog a Q(l + 



where f = - 2 Pw a ) 2 + 4|p i, 10 | 2 . Then, D and CC 
are determined by the minimum of the two expressions 
in fll|. 

For a — \, we find that S^(p A \{u B }) < <S ,(1 Hp.4|{n B }) 
is always satisfied and thus »S''- 2 HPA|{n B }) is the entropy 
used for computing the correlations. 

The concurrence for the density matrix (|17p takes the 
simple form C(p) = 2max{0, |poi,io|}- The entanglement 
of formation is then directly computed from Eq. (1111) . 

In Figure 2] we show the evolution of these three corre- 
lations as a function of the normalized time r = (r + 7)t. 
The exponential decay is due to the particular choice 
of maximally entangled initial states a = 1/2 for (a) 
the symmetric state ip = (main), and (b) the anti- 
symmetric state tp = 7r (inset). Although the qualitative 
behavior of the inset is similar to that of the main graph, 
we notice that such state evolution has a much slower de- 
cay rate of the correlations; this happens at the new rate 
T — 7. This result is in agreement with what would be 
expected for sub-radiant states and explains the higher 
value of the correlations at all times, when compared to 
the initial super-radiant ip = state. 

The above analysis is extended to all initial states 
< a < 1, see Figs. [5] and [6] We also find, for iden- 
tical qubits, that quantum correlations are stronger and 
more robust to decay than classical correlations. The dy- 
namical behavior is richer than the exponential decay of 
Fig. [4] because different initial conditions for the den- 
sity matrix allow for different roles for the corresponding 
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FIG. 6. (Color online) Quantum discord dynamics as a func- 
tion of the initial state preparations (a) a v= o, and (b) ol v — % 
states (main), a = 1/2 give the super-radiant and sub-radiant 
states, respectively. The inset notations are as in Fig. [S] for 
a = 0. 7 = 0.91 T, and V = 2.03 T. 



FIG. 5. (Color online) Quantum discord dynamics as a func- 
tion of the initial a v =o state preparations. The inset shows 
the dynamics of correlations D(pab) (solid blue), EoF(pab) 
(dashed, red), and CC(pab) (doubly dashed, green), for the 
(separable state) a — 0-plane. There is no collective damping: 
7 = 0, and V = 2T. 



collective parameters. For a ^ 1/2, the last terms of ma- 
trix elements poi,oi an d Pio.icb an d the imaginary part 
°f Poi.iOj are non-zero and hence the dynamics exhibits 
oscillations depending on the value of the d-d interaction, 
which is weak for the case of the diluted molecules first 
considered; accordingly, we consider in the first place the 
scenario of no collective decay due to 7 (decay contri- 
butions are only those giving by T), as shown in Fig. [5] 
for the initial state dynamics 01^=0. We also consider, 
in Fig. |51 the situation where 7 7^ plays an important 
role, for the initial state preparations (a) a v= o, and (b) 
ct v= Tr- The correlation dynamics for these two cases is 
shown in Figs. [5]and[f)] 

When the collective decay rate is not effectively present 
in the dynamics, as, for example, in the case of diluted 



molecules in a dispersive environment for which their col- 
lective rate is much less than both the spontaneous indi- 
vidual emission rates and the dipolar interaction [l9|, [22| ; 
or in a pair of emitters that interact with a plasmonic 
channel where the collective stregths V and 7 can effec- 
tively be switched on and off as the distance between the 
emitters is varied, a mechanism that produces an effec- 
tive phase difference between V and 7 [29| . 

The main graph in Fig. [5] plots the quantum discord 
dynamics as a function of a for two different set of initial 
state preparations; the same qualitative behavior is ob- 
served for CC and EoF (not shown). In Fig. [S]we see that 
oscillations due to the d-d interaction begin to appear 
and to coherently affect the dynamics of the quantum 
discord for a values closer to or 1; this is so because in 
these cases the amplitude of the trigonometric functions 
in the matrix elements Eq. (|16p are closer to their max- 
ima. The inset in Fig. [5] shows this for the discord, the 
EoF and CCs for the plane a = (separable state 1 10} ) . 
For the maximally entangled state (a = 1/2) all the cor- 
relations have a monotonic exponential decay, indicating 
an independence from the interaction V. For other initial 
states set by a, this decay exhibits oscillations which in- 
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crease their amplitude as a tends to the separable states 
1 01) or 1 10). The oscillations eventually vanish and the 
correlations decay to for all a. Although the main 
graphs of Fig. [5] only plot the discord, our calculations 
show that the correlations hierarchy CC < EoF < QD 
holds for all time t and for all a (see the inset for the 
case a = 0). 

The case 7^0 (Fig. [5]) plays an important role in 
the system's evolution. Here, two emitters interact via 
a dipole force throughout all modes of the electromag- 
netic field in the vacuum state: the collective (incoher- 
ent) decay 7 (|15p has to be taken into account, and we 
introduce the full expressions given for 7 and V, and 
their relation to the Einstein coefficients, dipole orienta- 
tion, and the dispersive medium, as given by Eqs. (fT4]) 
and (fT5)) . When compared to Fig. [5J the effect of the 
collective decay rate translates into a much longer time 
for the survival of correlations. In fact, the inset, which 
is common to Figs. |6ja) and (b), shows that quantum 
correlations decay very slowly and thus last for a much 
longer time than that obtained in the absence of collec- 
tive decay in Fig. [5] This said, in Fig. [6] there is a 
clear distinction between the initial (a) a v= o, and (b) 
aip=Tr states. In the latter, the correlations last for a far 
longer time before their final decay; in particular, the 
case of the anti-symmetric a v=T = 1/2 state exhibits a 
very slow decay of the discord, as compared to the case of 
the symmetric a v= o = 1/2 state (note the difference in 
the time scale of the two graphs). This is also reflected in 
the inset of Fig. EJb), for the case a = 0. The difference 
in physical behavior between graphs (a) and (b) owes its 
origin to the nature of the super-radiance phenomenon, 
which is manifest here in graph (b) for a = 1/2. In any 
case, in both graphs, we have quantum correlations al- 
ways greater than classical correlations. 

By computing the correlation dynamics for two emit- 
ters of identical dipole moments (parallel to each other), 
and perpendicular to the interqubit distance i"i2, basic 
features can be identified. The relevant values of the in- 
teraction energy and collective decay rate are given by 
Eqs. (rfj} and (|T5j) : in Fig. |we consider 7 = 0.91 T 
and V = 2.03 r, which corresponds to a distance 772 = 
0.108 Ao- Even though we now have less oscillations in 
the discord, a higher degree of correlations than those 
of Fig. [5] is obtained due to the presence of collective 
damping. Note that for initial configurations of the type 
a —> 1 or a — > 0, quantum discord and entanglement hold 
a greater value than classical correlations and, as can be 
seen from Figs. [5] and [6l the discord is the correlation 
that most benefits due to this collective effect. The en- 
tropy condition Eq. (fT9")) is always satisfied, and hence 
quantum correlations continue to be always greater than 
classical correlations: CC < EoF < QD, for all a. 

According to the above inspection, a calculation about 
the way the collective parameters influence the dynamical 
behavior of correlations is in order. From Eqs. (fT5|) it is 
clear that a variation of V produces a change in the oscil- 
lations of the density matrix elements and therefore in the 




FIG. 7. (Color online) Correlations D, EoF, and CC, as 
a function of time and the interqubit distance, the latter 
normalised by the emitters characteristic wavelength Ao = 
27r/fco = 2-kc/uq. The initial state preparation corresponds 
to a = 0. 



correlations dynamics: the degree of correlations is in- 
creased as the collective damping is augmented. Figure [7] 
shows the behavior of the three relevant correlations as a 
function of the normalized distance 772/ Ao, in the regime 
where V goes from 0.2r to 2.6r, and 0.1 < ^ < 0.9, or, 
equivalently, for 0.1 < y 2 - < 0.4. In Fig. [7] we consider 
the initial condition a = 0, following the results shown in 
Figs. [S]and|ni as such a state makes the strongest case for 
the influence of the collective parameters. As expected, 
a signature in the correlations dynamics shows up after 
the interaction strength becomes either of the order of or 
greater than the decay rate T, which occurs for ^2. < 0.2. 
As in the previous case, the most robust correlation is 
the quantum discord and the classical correlation is the 
least assisted by the collective effects throughout the en- 
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FIG. 8. (Color online) Dynamics of correlations for emitters 
under resonant laser field excitation: vl = Vo = v\ ii = £■ 
Notations for D, EoF, and CC are as in previous figures, (a) 

I = 1.5 r, 7 = o.9i r, v = 2.03 r ; (b) i = 10.0 r, 7 = 0.91 r, 

V = 2.03 T; (c) i = 10.0 I\ 7 = 0.97 T, V = 10.45 F. Fhe 
initial state is the super-radiant a v= o = 1/2. 




n 

FIG. 9. (Color online) Emitters correlations (notation as in 
Fig. [5| under resonant laser excitation vl = v. Parameters 
are as in Fig. [8] Initial state: antisymmetric sub-radiant 
av=ir = 1/2. 



VI. EMITTERS CORRELATIONS UNDER 
COHERENT LASER EXCITATION 



tire dynamics. Hence, the degree of correlations can be 
dynamically enhanced when the collective (incoherent) 
damping is of the order of the (weak) interaction energy 
and comparable to the spontaneous decay rates of the 
single quantum emitters. 

Lindblad [IH conjectured that for a general density 
matrix p, classical correlations are greater than half of the 
mutual information and as such, they are always greater 
than quantum correlations. He showed that for an inco- 
herent mixture state, CC(p) = I(p), and that for a pure 
state, CC(p) = \l{p) 1 it seemed natural to conjecture 
that CC(p) > \l(p) Hf|. Simple numerical examples, 
for several different configurations of the constants hi in 
Pm , confirm that for the cases of incoherent mixtures and 
pure states, this result holds. There are states for which 
the conjecture does not hold, however: it is straightfor- 
ward to show that for the set of Bell states \hi\ = 1 (pure 
correlated states), I{pBeii) = 2, hence CC(pBeii) = 1- 
On the other hand, for hi = hi = 0, and \h 3 \ < 1 (inco- 
herent mixture), CC(p lncoh ) = I(pincoh)- For example, 
for h 3 = 0.6, CC(p incoh ) = I(pincoh) = 0.278. Finally, 
the choice hi = h 2 = 0.8, h 3 = —0.6 gives I(pm) = 1-078, 
and quantum correlations D(pm) = 0.547 > CC(pm) = 
0.531. 



We now turn to the external all-optical control of the 
quantum emitters correlation dynamics. We consider the 
case of a continuous laser field ^coswit locally applied 
to each of the emitters, where ii denotes the amplitude 
of the laser field acting on qubit i, and ojl is the laser 
frequency. Without loss of generality, we consider qubits 
of the same transition frequency u>, under the resonant 
condition ujl = luq = uj. 

The solution of the master equation describing the dy- 
namical system is fully numerical because the structure 
of the density matrix, when taking into account the laser 
field, does not admit a general analytical solution. 

Figure [5] shows the results for initial state preparation 
in the super-radiant ^= (| 01) + 1 10} ) state. For external 
laser field amplitude below the interaction energy V (Fig. 
IHJa), and 7, and V as in Fig. [5]), the emitters dynamics 
exhibits a change in the correlations hierarchy and, after 
a brief transient period, quantum correlations go below 
the classical ones. This change in physical behavior, for 
weak laser excitation, occurs in almost the whole dynam- 
ics. Even though such laser excitation can overdamp the 
evolution of the quantum correlations well below the clas- 
sical correlations, such evolution can be controlled and 
the degree of the quantum correlations can be optimised 
by appropriately tailoring both the amplitude of the field 
and the dipolar interaction between the qubits; the latter 
being obtained by varying the interqubit distance. The 
results plotted in Fig. IHa), for the initial state prcpara- 
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tion in the sub-radiant -1=(|01) — 1 1 0) ) state arc in clear 

contrast with those of Fig. HJa), and, although quantum 
correlations are no longer greater than classical corre- 
lations, they become of the same order throughout the 
entire dynamics. The early stage disentanglement (ESD) 
experienced by the emitters in Fig. M,&) is delayed for 
about ten times that of Fig. IBJa)- It is worth pointing 
out that even though there is no entanglement after ESD 
occurs, from this point on there arc non-zero classical 
and quantum correlations present in the evolution. The 
order hierarchy for the anti-symmetric state is such that 
QD > CC at all times. This is not always the case for 
the symmetric state, and the hierarchy depends on the 
interplay between V and I. 

In Figs. [5Jb) andJSJb) we have kept the interaction en- 
ergy constant (fixed collective parameters V, and 7) and 
increased the value for the laser amplitude well above 

V (I = 10r). Now the discord evolves as the most ro- 
bust correlation. For the case of the symmetric state, 
however, not always QD > CC (initially, CC > QD), 
whereas for the anti-symmetric state, we always have 
QD > CC. The dynamics becomes more complex with 
the inclusion of the laser field and a dependence with the 
d-d interaction now appears, in contrast to the case of 
no laser excitation. This is seen in Figs. HJb) and (c): 
for a fixed laser amplitude, the variation in the d-d inter- 
action produces a qualitative change which favours the 
degree of discord, above classical correlations and entan- 
glement, as V increases. The early stage disentanglement 
doesn't show major changes. Now we turn to Figs. IHJb) 
and (c): quantum correlations are always higher than 
classical correlations and the effect of an increment in 

V noticeably prolongucs the lifetime and the degree of 
correlations. For example, ESD in graphic) is around 
50 times that of [H^c) (notice the difference in the corre- 
sponding time scales). 

In short, the same laser field can reverse the behavior 
observed in Figs. [SJa) and EJa) by appropriately chang- 
ing the interqubit distance and by shining the laser with 
higher intensity. This control procedure replicates the 
earlier result of quantum correlations being above classi- 
cal correlations. It is also clear that at all times, quantum 
correlations are far greater than entanglement. The pres- 
ence of the external field perturbs the asymptotic decay 
of entanglement and in turn produces an early stage dis- 
entanglement phenomenon |2lL |27| which is verified in all 
the graphs of Figs. [8] and [9] 

For completeness, in Fig. [TU]we compare the graphs for 
a separable state as the initial condition, and for (a) 'ab- 
sent' 7 and (b) 'active' 7. The correlations dynamics has 
been plotted for V = 10.45 T, with laser amplitude I ~ V, 
following the parameters of graphic) which give a high 
degree of discord for symmetric states. Unlike in the pre- 
vious cases, we now obtain oscillations in the discord and 
entanglement. Stationary correlations, which are null in 
the absence of collective damping (see Fig. ITUT a)). can 
now be obtained due to the presence of damping 7 in 
Fig. ITUl' bh Interestingly, as the correlations stabilise, a 




FIG. 10. (Color online) Correlations dynamics under the ef- 
fect of an on-resonance external field uil = cjo, with amplitude 
£ = 10 T; V = 10.45 T. (a) Emitters with collective damping 
'switched off'; (b) Emitters with identical and parallel dipoles, 
with 7 = 0.97 T. Initial state: a = 0. 

quantum discord more robust than classical correlations 
emerges. This kind of quantum information or correla- 
tions is of an entirely different nature to entanglement, 
which decreases to zero much earlier than QD and CC. 
We conclude that also for quantum emitters under laser 
pumping, quantum correlations can become greater than 
classical correlations and that this feature is more 'natu- 
ral' in antisymmetric sub-radiant states, thus disproving 
Linblad's conjecture [251 ], 

VII. SUMMARY 

We have given a thorough analysis of the dissipative 
dynamics of total, quantum, and classical correlations 
for a set of interacting quantum emitters characterized 
by their coupling energy and their collective and indi- 
vidual decay rates, and analytically demonstrated the 
conditions under which quantum correlations are always 
greater than classical correlations, and the role played 
by entanglement during such dissipative quantum evo- 
lution. We have shown that it is the entanglement of 
formation the metric to be used as the emmiters' quanti- 
fier of entanglement and that the use of concurrence for 
this purpose can give misleading results about the actual 
degree of generated entanglement. We have also shown 
that the dynamics of correlations followed by the emit- 
ters can be controlled and engineered in the presence of 
a coherent laser field. We have quantified and illustrated 
the interplay between entanglement and quantum discord 
stressing the importance of the persistent quantum corre- 
lations in several different scenarios where entanglement 
vanishes. 

In cases where classical correlations are greater than 
quantum correlations, a coherent control mechanism can 
be implemented so that the latter are enhanced and be- 
come greater than the former. This result is corroborated 
for driven emitters in the case of strong laser-dipole cou- 
pling: by preparing the initial symmetric state, once en- 
tanglement vanishes there is a rapid stabilisation of corre- 
lations and the discord becomes greater than the classical 
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correlations. Moreover, for strong laser excitation, now 
starting with the antisymmetric state, the discord is, at 
all times, greater than the classical correlations, regard- 
less of the emitters' entanglement content. This result is 
resilient to scenarios where the collective (incoherent) de- 
cay rate reflects a strong coupling of the emitters to the 
surrounding environment. In fact, the strength of the 
quantum correlations can profit from large collective de- 
cay rates, if the emitters strength coupling and the laser 
field excitation are appropriately tuned. 

From an experimental viewpoint, set-ups comprising 
sing le molecules 0, [22| or semiconductor quantum dots 
[261 ] could be employed to demonstrate our findings. In 



particular, single-molecule techniques that allow fem- 
tosecond pulse-shaping are currently available, and con- 
trolled coherent superpositions and basic femtosecond 
single-qubit operations in single organic molecules have 
been produced and manipulated even at room tempera- 
ture (22( 1 - 
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